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Abstract 



We study the envy- free cake-cutting problem for d + 1 players with 
d cuts, for both the oracle function model and the polynomial time 
function model. For the former, we derive a OK^)^ 1 ) time match- 
ing bound for the query complexity of d + 1 player cake cutting with 
Lipschitz utilities for any d > 1. When the utility functions are given 
by a polynomial time algorithm, we prove the problem to be PPAD- 
complete. 

For measurable utility functions, we find a fully polynomial-time 
algorithm for finding an approximate envy-free allocation of a cake 
among three people using two cuts. 

1 Introduction 

Suppose you have a cake represented by the interval ([0, 1]), and you would 
like to divide it among n persons fairly. Each person may have a different 
opinion as to which part is more valuable. There is a big literature on this 
problem in economics, political science and computer science |24[ [25l [T2[ [31 
13 [20], [19] [H [30]. In particular, it is proved using a fixed-point argument 
that this problem has an envy-free solution [26] [28| |2~7] . In other words, it is 
possible to cut a cake into re pieces {X = {Iq, ii, • • ■ , l n -i} from left to right 
along [0, 1]) using n — 1 cuts and to allocate one piece to each person (player 
i assigned piece l^u\ for permutation n(i)) so that everyone values his or 
her assigned piece no less than any other piece. The question is: is there an 
efficient algorithm that finds such a cut (called (re — l)-cut subsequently) of 
the cake? 

A related but less demanding solution than an envy-free solution is that 
of proportional cuts. That is, each person gets a piece which he or she values 
more than \jn of total. Its complexity has recently completely solved by 
Edmonds and Pruhs |10j|llj . 

For the d + 1-person envy-free cut problem with exactly d cuts under 
our consideration, however, progress in complexity analysis has been limited. 
The existence of such a solution was proven by Stromquist [26] with a fixed- 
point argument. His proof implies that an e-approximation can be found 
in time exponential in input size O(log^). Let N = - throughout our 
discussion. 

We establish three main results: 1) When the best choices of the play- 
ers are given by polynomial-time algorithms, we prove that the problem is 
PPAD-complete. 2) If the choices are given by a functional oracle, we derive 
a ) matching bound for the query complexity of cutting a cake for 
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d + 1 players. Despite a strong connection, mathematical and complexity- 
wise, between equilibrium computation and fixed point computation, this 
is the first matching query complexity result for an equilibrium computa- 
tion problem. We know of no such results for Nash equilibrium, which is 
also PPAD-complete with a strong tie with fixed-point computation. 3) For 
the special case of measurable utility functions, we make a simple observa- 
tion: there is a fully polynomial-time approximation scheme for finding an 
approximate envy-free allocation 2-cut of a cake among three people. We 
sketch our approaches as follows. 

PPAD completeness: First, we capture the concept of approximation 
by defining a discrete envy- free cut (set) (tt, l(°),lW, • • • ,XW) such that 
player i prefers the 7r(i)-th piece of the d-cut X^) for some j, and the d- 
cuts {X(°\ X^, ■ ■ ■ ,1^} are within e distance (in metric) of each 
other. Such a solution converges to an exact envy- free cake cut as the 
distance bound for the d-cuts goes to zero. Using barycentric coordinates, 
the o?-cut can also be represented as (xo,x±, ■ ■ ■ ,Xd) with Iq = xo/N,li = 
xi/N, ■ ■ ■ ,ld = X d/X. Note that (xo, xi, ■ ■ ■ , x^) is a point on the standard 
c?-dimensional simplex. To prove that the problem is in PPAD, we reduce it 
to the problem of finding a fully colored base cell in a triangulated Sperner 
coloring of a cZ-simplex by using Kuhn's triangulation |18j . This can be done 
by a two-stage process: labeling and coloring. First, for a (i-simplex and a 
Kuhn's triangulation with vertex set V, a labeling C : V — > {0, 1, ■ • • , d} is 
valid if \fX,Y e V and X,Y on the same base cell, C{X) ^ C(Y). Then 
for any labeled vertex X, we define a coloring C : V — > {0, 1, • • • ,d} such 
that C(X) = i if player C(X) prefers the i-th piece of the cut X. By a mild 
condition for the utility functions, C is a proper Sperner coloring. Therefore, 
the key point here is to find a polynomial time labeling rule. We define the 
labeling rule as: C{X) = Y2i=o m °d (d + 1) with a proof of its validity 
for Kuhn's triangulation. On the other hand, we design a reduction based 
on the 2D BROUWER problem [2H El U\ [8] for its proof of PPAD-hardness. 

Matching bound in the Oracle function model: We derive a #((7) d_1 ) 
time matching bound for the query complexity of cake cutting for d + 1 
players with Lipschitz utilities. The tight upper bound requires a divide- 
&-conquer method that finds a balanced cut of the simplex. It is made 
possible by Kuhn's triangulation of the simplex, and our labeling method 
for the envy-free cake cutting problem that allows an efficient parity checking 
of the boundary. For the lower bound, the results are obtained by a reduc- 
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tion to the zero-point problem \14\ [T5l [161 [6]. The reduction is achieved in 
two steps. First, we reduce the zero-point problem for direction-preserving 
functions to the problem of finding a discrete fixed-point on a hypergrid. 
In the second step, we prove that the hypergrid can be embedded into the 
original d-simplex for cake-cutting such that its coloring can be extended to 
a proper Sperner coloring of vertices in the triangulated simplex. 

Instrumental to our matching bound for envy-free cake-cutting, we prove 
a matching bound for the SPERNER problem for any constant dimension 
d > 2, in the oracle function model. This was an open problem, while for 
the case of d = 2, a tight bound was known by a lower bound of Crescenzi 
and Silvestri [9] and an upper bound of Friedl, et al [13] . This matching 
bound for the SPERNER problem may have other applications for fixed- 
point based solutions. 

Fully PTAS for three players with measurable utility functions: 

Finally, for the special case of measurable utility functions, we are able to 
utilize their monotone properties to construct an e-approximate envy-free 
solution in time polynomial in log(^) for three players. 

We still rely on the general approach of branch-&-bound on parity but 
exploit the monotonicity of the best choice along certain lines of the possible 
cuts to make an efficient count of the index along the boundary. Any player 
with a measurable utility function would prefer A to B for B C A. Therefore, 
when one cut is fixed, one of the three pieces is fixed. Any player's preference 
on the other two pieces will change monotonically as another cut changes 
from left to right. Using the barycentric coordinate X = (xo,xi,X2), along 
the line of fixed xo, let x\ increases from to N — xq. The preference of a 
player's choice will start with the last piece I2, to lo, and to l\ (one or two of 
them may be missing). Similar monotone property holds when X2 is fixed. 
For each player, we can find the boundary point along those lines by binary 
search and the break point of the choices will can be obtained in 0(log N) 
time. 

We cut the space along those two directions, so that the choice function 
of each individual player will be monotone along those directions. Because 
of the monotonicity and using the above procedure, we can calculate the 
indices of edges along those lines efficiently. Therefore, the indices of the 
two regions split by the cut will be decided quickly. We will stay on the 
region with an odd index so that we should end with one that is a diamond 
shape polygon consisting of at most two base cells. Because of parity, one of 
two cells is a fully colored base triangle. The overall query complexity and 
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running time will be of 0(log 2 N). 



2 Triangulation and Index 

Our results are based on Sperner's Lemma and its generalizations using the 
concept of the index of a region [29]. These results have been fundamental in 
discrete fixed point computation (see, e.g., [22] and [29] ) and establishing 
Brouwer's fixed point theorem [JJ. 

Starting at two dimensions, a triangular grid of scale 1 is an ordinary 
triangle A which has three vertices and one base cell. A triangular grid of 
scale N places N-l equally spaced line segments parallel to each of the three 
edges of A and divides the triangle into iV 2 base cells. 

We refer to the three vertices of A as corner vertices, and denote them 
by Do, Di, and Di- The vertices along the edges of A are referred to as 
boundary vertices. Other vertices are referred to as internal vertices. Edges 
of a base cell are referred to as base edges. Each vertex x = (i x Dq + j x 
D\ + k x D2)/N is represented by (i,j, k) where i, j, k > 0, i + j + k = N. 
We call it the barycentric coordinates of the vertex. 

By barycentric coordinates, Do is represented by (N, 0, 0); D\ by (0, N, 0) 
and Z?2 by (0,0, N). Boundary vertices along Do and D\ are the ones in 
the form (i,j, 0) with i,j > Q,i + j = N. Other boundary vertices are 
defined similarly. For any interior point represented by (i,j,k) we have 
i,j,k>0,i+j + k = N. Let V = {(i,j,k) : i, j, k > 0,i + j + k = N}. 

Base cells in the triangulation are oriented in the clockwise order of their 
vertices and base edges are oriented according to the clockwise order of their 
base cells. See Figure [T] in Appendix A.l. 

A coloring (j) : V — > {0, 1, 2} is a Sperner coloring if and only if for any 
vertex x = (xo, xi, X2), 0(xo> £2) = j € {0,1,2} implies Xj > 0. Sperner 
Lemma states that a triangulated triangle with a valid Sperner color has a 
base cell such that its three vertices have different colors. 

Given a Sperner coloring <f> : V — > {0,1,2} of all vertices in V, let 
sign(5, (j)) and sign(e, 5, cj>) denote the sign of a base cell 5 and the sign of 
a base edge e in 5 respectively. The sign of e = (u, v) is 1 (or -1) if the 
colors of its two vertices are and 1 and the orientation of e in 5 is from 
color vertex to color 1 vertex (or from 1 to 0). We denote the sign of 
a base edge by sign(e, <ft) if there is no ambiguity on its orientation. In 
all other cases, sign(e, (j)) = 0. The sign of a base cell 5 is defined to 
be the sum of the signs of its three base edges. Therefore, sign(5, eft) = 
sign(e\,5, <p) + sign(e<2,5, (ft) + sign(e%, 5, (/>). 
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We may verify the following by a simple case analysis. 

Proposition 1. For any base cell, its sign is 1 (or —1) if and only if its three 
vertices are colored with 0,1,2 in the clockwise (counterclockwise) order. In 
all other cases, its sign is zero. 

The index of a connected set of base cells A, with respect to the color <p 
is defined as: 

index(A, 0) = ' s ^{sign{5, <f>) : 5 a base triangle G A} 

Lemma 1. [29] For a triangulated triangle A with colors (ft : V — > {0, 1, 2} ; 
there are at least \index(A, <j))\ fully colored base cells. The index can be 
calculated by summing the signs on its boundary base edges. 

See Appendix A. 2 for the proof. The result also holds for general poly- 
gons in 2D. 

To generalize the same result to a higher dimensional polyhedron P, we 
consider a simpler version of index that is defined mod 2 and is used in [8]. 
For a d-dimensional simplex with vertices assigned d + 1 different colors 
{0, 1, • • • , d}, we define its index as 1. Otherwise, it is defined to be zero. Let 
V(P) be the vertices of its triangulation. With respect to a color <p ■ V(P) — > 
{0, 1, • • • , d}, its index is defined as index(P, <p) = ^2g eP index(5, <f>), where 
(5's are e?-dimensional simplices in the triangulation of P into simplices. 

Denote by dP the boundaries of P. Note that the triangulation of P 
induces a triangulation of dP into (d — l)-dimensional simplices. We define 
index d -i(dP,cp) = Y^8 d _ ie dP index d-i(^d-i,(p)- 

We need the following discrete version [6] of standard results on the index 
defined here. 

Proposition 2. index(P,(f>) = index d-i(dP,(p) mod 2. 
See Appendix A. 3 for the proof. 

2.1 Kuhn's Triangulation 

In this section, we briefly introduce Kuhn's triangulation [18] for a simplex. 
Kuhn's triangulation has the advantage of being a balanced triangulation 
and it helps us derive a much improved algorithm for the envy-free cake 
cutting problem. 

Let us start by explaining the Kuhn's triangulation of a unit cube in 
d dimensions. Let vq = (0, 0, • • • , 0)i X d be one of the corners of the cube. 
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The diagonal vertex to it would be Vd+i = (1, 1, • • • , l)ix<2- Suppose e, is a 
(i-dimensional unit vector such that en = 1 and eij = for all i ^ j. Kuhn's 
method partitions the cube into d\ simplices. Let tt := (vr(l), vr(2), ■ • • , 7r(d)) 
be any permutation of the integers 0, l,--- , d — 1. Each permutation tt 
corresponds to one small simplex whose vertices are given by v\ = 
v%r x + e n{i) and = v . 

These simplices all have disjoint interiors and their union is the d-cube. 
It is not difficult to verify this, since any vertex x = (xo,x%,- ■ ■ ,Xd-i) is 
an interior point of A^ if and only if 1 > a^m > x^m) > • • • > ^(d) > 0. 
Appendix A. 4 illustrates Kuhn's triangulation on a 3-cube. 

Now, we are ready to explain the Kuhn's triangulation of a simplex. Let 
N be an integer bigger than 1. Take a unit d-cube and use parallel cuts 
of equal distance to partition it into N d smaller d-cubes of side length 
Then, partition each small cube into d\ simplices using the above method. 
Now, observe that the unit cube can also be partitioned into dl big simplices 
first and each big simplex contains N d smaller simplices or base cells. The 
proof of the consistency of the two processes can be found in Appendix A. 5. 

Based on the equivalency of the two partitioning processes, we choose one 
of the big simplices, for example the one corresponding to tt = (0, 1, 2, • • ■ , d— 
1) and the smaller simplices that are contained in that will define its trian- 
gulation. A vertex X in this big simplices can be represented by barycentric 
coordinates X = (xo,xi, • • • ,Xd) by a transformation as illustrated in Ap- 
pendix A. 6. More details of the transformation can be found in page 42 of 

We have the following property of the triangulation: 

Lemma 2. For given X = (xo, x±, ■ ■ ■ , Xd) and Y = (yo, yx, ■ ■ ■ , yd), define 
5x-Y = m axvig{o,i,--- ,d}{\ x i ~ yi\\- If A Y are in the same base cell in 
Kuhn's triangulation, then Sx-Y = 1- 

See Appendix A. 7 for the proof. 

3 Finding a Sperner Simplex under Oracle Func- 
tion Model 

In the oracle function model, the function value at a point (or color of the 
point) is given only when it is queried and it remains the same when further 
queries are performed on the same point. 

We prove that the oracle complexity of finding a Sperner base simplex 
under oracle function model in d dimensions is of 9(N d ~ l ). Such a matching 
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bound was known only for finding a Sperner's fully colored base cell in a 
two dimensional N x N grid [9] |13j . Our extensions into higher dimensions 
make use of the methodologies originated in the zero point computation on 
hypercubes. 

The matching bound derived in this section will be essential for solving 
the envy-free cake-cutting problem. It could also be of independent interest 
for other equilibrium problems that are based on fixed point computation. 

To derive the upper bound, we define the concept of balanced triangu- 
lations: 

Definition 1. A simplex P is triangulated into balanced simplices of gran- 
ularity g = if 

1. P is fully contained in the unit cube [0, l] d ; 

2. every parallel plane along the coordinates X{ = g x j (i = 1, 2, • • • , d, 
j = 0, 1, • • • , N) cuts through P along the facets of the base cells of the 
triangulation, i.e., the parallel plane will not cut into the base cells; 

3. the number of the d- dimensional simplices of the triangulation within 
any cube of side length g is constant. 

By construction, Kuhn's triangulation is a balanced triangulation. 

Lemma 3. For any balanced triangulation, there is an algorithm that finds 
a Sperner base cell in time 0(N d ^ 1 ) if vertices are colored by a Sperner 
coloring. 

Proof. We fit the balanced triangulated simplex P into the unit cube [0, l] d 
which guaranteed by condition 1 of Definition 1. By condition 2, we use 
parallel plane along the coordinates to cut the cube. 

For the correctness, we note that, which sub-hypercube to look into will 
be determined by applying Proposition [2] on the triangulated simplex. As 
the boundary conditions give an odd index for the initial simplex because 
of the valid coloring, each time one of the two parts of the cut simplex will 
be odd. The procedure can proceed until the last base cube. Then, we can 
simply examine up to C remaining simplices contained in this base cube, 
where C is a constant (note that because of d is a constant, the function: 
f(d) is a constant function). 

For complexity, our algorithm will be doing a binary cut along the d 
coordinate one after another. It takes d such cuts to reduce a hypercube 
into half of its size(in length). Therefore, in d\og 2 (|) cuts, we reduce the 
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unit hypercube into a base cube of side length g. We upper bound the 
time complexity with the total time necessary for the hypercube (the actual 
number of operations will be less on the simplex and its triangulation) . As 
the size reduces geometrically, the total number of operations is dominated 
by the number of operations we do at the first d cuts. For each cut, we need 
to apply Proposition [2] to calculate the index of the boundary simplices of 
(d — l)-dimension, which requires a computational time 0((|) rf_1 )- The 
computational upper bound follows. □ 

We also derive a similar lower bound: 

Lemma 4. For any algorithm that finds a Sperner base cell for any tri- 
angulation of a d- dimensional simplex with a Sperner coloring, there exists 
some input triangulation such that the algorithm takes time Q(N d ~ 1 ). 

The proof of the above lemma is build on a deep result of Chen and 
Deng [6]. It is explained in details in Appendix A. 8. 

Lemma [3] and Lemma [H result in a matching bound as follows. 

Theorem 1. (matching bound) Given a balanced triangulation where all 
vertices are colored by {0, l,--- ,d} by a Sperner coloring, a Sperner base 
cell can be found in time 8(N d ^ 1 ). 

4 Envy-Free Cake Cutting and Sperner Lemma 

As far as we know, the first proof of the existence of envy-free cake cutting 
solutions using Sperner 's lemma is by Simmons [28]. Su [27] uses a similar 
argument to develop a computational procedure to derive an approximate 
envy-free cake cutting solution, by a labeling process on barycentric sub- 
divisions of a simplex (See figure [Has in Appendix A. 9 as well as |2j[27j). 
However, Su's method creates simplices with large aspect-ratios that make 
the process converge rather slowly. Instead, we use Kuhn's triangulation. 

4.1 Utility functions and envy- free solutions 

Consider a set / = {0, 1, • • • , d} of d + 1 players. Each player i E / has a 
utility function Ui defined on the Borel space of the line segment L = [0, 1]. 
Our utility functions are required to satisfy the following two conditions: 

• Nonnegativity condition: Uj(0) = and Uj(^ 0) > 0. 

• Lipschitz condition: For any interval [x,y] C L, Ui([x,y]) < K x \y — x\. 
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We use d cuts to partition L into a set S of d + 1 disjoint segments of 
lengths Iq,1\, ■ • • , such that Yli=o = Using barycentric coordinates, 
we restrict our discussion to integer vectors (xq,x±, • • • ,Xd) such that Iq = 
xq/N, l\ = xi/N, • • • ,ld = xj/N, with Yli=o Xi = ^- All possible partitions 
form a d-dimensional simplex A d with d+1 vertices. The i-th vertex of A d 
is represented as Nei, where i = {0, 1, ■ ■ ■ , d} and is the unit vector whose 
i-th coordinate is 1. 

By the nonnegativity condition of utility functions, every player will 
strictly prefer the nonzero segments to the zero segments. Hence, we have 
the following boundary preference condition. 

Property 1. Boundary Preference Property: consider any boundary vertex 
X that belongs to a boundary which is incident to the i-th corner but not the 
j-th corner. In the cake cutting defined by X, every player strictly prefers 
the i-th segment to the j-th. 

The above property actually ensures a Sperner coloring. We can now 
give a sketch of the envy-free cake cutting problem using this connection. 
The argument has two stages: labeling and coloring. 

For the simplicity of exposition, consider the problem for three players. 
For the case of 3 players, the closed set of all possible cuts is a triangle. 
As in the previous section we place N — 1 equally spaced line segments 
parallel to each of the three edges of the triangle and divide it into iV 2 
equal-size base triangles. Let V be the set of vertices of all base triangles, 
i.e., V = {(x , x\,x 2 ) : xq, X\,x 2 > 0, xq + x\ + x 2 = N}. 

Next, we partition V into three control subsets Vq,Vi,V 2 . Starting by 
assigning (N,0,0) to V , {N - 1,1,0) to V ± , and (N - 2,2,0) to V 2 . The 
rest of V is partitioned in such a way that the three vertices of each base 
triangle belong to different subsets Vt's. This can be done by defining V t = 
{(xo, Xi, x 2 ) ■ x\ — x 2 = t(mod 3) for xq, x±,x 2 > 0, xq + x\ + x 2 = N}. 

Now, we should color V. For any vertex (xq, x\, x 2 ) G Vt, we let player 
t choose, among three segments, [0, Iq], [Iq, Iq + l±], [Iq + l\, 1], of /, one that 
maximizes his utility. For simplicity of presentation, we should assume a 
non-degenerate condition that the choice is unique. The general case can be 
handled with a careful tie-breaking rule. If the optimal segment is the one 
of length l s , < s < 2, we assign color s to the vertex (xq, x±,x 2 ). We claim 
that the above coloring is a valid Sperner coloring. This can be easily checked 
by the assumption of utility functions. Since Uj(0) = 0, the three vertices 
(N, 0, 0), (0, N, 0), (0, 0, N) of the large triangle must be colored by 0,1,2 
respectively and the vertices on the edge (N, 0,0) — > (0,iV, 0), (0,iV, 0) — > 
(0,0, AT), and (0,0, N) -> (N,0,0) will be colored either by or 1, 1 or 2 
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and 2 or respectively by our coloring procedure. Hence, it satisfies the 
boundary condition of Sperner lemma, and the coloring is valid. See Figure 
[5] in Appendix A. 10 for an example. 

Since the three vertices of each base triangle belong to three different 
subsets, or we say three different players, if we find a fully colored base 
triangle, then on the three vertices of this triangle, different players prefer 
different segments. By Sperner Lemma, there exists at least one fully colored 
base triangle. By refining the triangulations, the fully colored base triangles 
become smaller and smaller, and a subsequence of the base triangles will 
converges to a fixed point. Such a fixed point is an envy free solution for the 
cake cutting problem. Therefore, there always exist an envy-free solution 
for 3 players case. 

The case of d + 1 players for d > 2 can be handled in a similar way. 
In that case, the closed set of all possible partitions of the cake forms a 
(i-dimensional simplex. 

In Kuhn's triangulation the labeling can be done by using barycentric 
coordinates. X = (xq, x%, ■ ■ ■ , x^)- For any vertex X of the base simplex, let 
W(X) = Y^i =0 m. We assign vertex X to subset V t if W(X) = t(mod d+l). 
This will partition the vertices of the base simplices into d+l control subsets 
Vo,Vi, - ■ ■ ,Vd- This is a suitable labeling because 

W(vi) = Wivir 1 ) + ((tt(0 + 1) - Tr(i)) = wtytr 1 ) + I- 

This also proves the existence of a d-cut solution. 

Theorem 2. (26), There is an envy-free cake cutting solution for d+l 
players that uses only d cuts. 

It is not hard to see that a fully colored base simplex represents an 
approximate envy-free cake cutting solution. In fact, by using the Lipschitz 
condition defined above, one can find a cake cutting solution with maximum 
envy e through a triangulation in which the sizes of all base simplices is 
bounded by e/K. Motivated by this observation, we define the discrete cake 
cutting problem and derive its computation and oracle complexity. 

5 Complexity of Discrete Cake Cutting 

To formalize the analysis for the cake cutting problem, we introduce a dis- 
crete version of the envy-free allocation of a cake among d + l people using 
d cuts (the d-cut problem for short). We use the barycentric coordinates 
(xq,x\,--- ,Xd) as in the previous section restricting x to integer vectors 
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satisfying Yli=o x i = N. We define two d-cuts x and y to be adjacent to 
each other, if Vi € {1, 2, • • ■ , d} : \xi — yi\ < 1. We call them afhne adjacent 
to each other if they are adjacent to each other and \xq — yo\ < 1. 

A discrete cake cut is defined to be a set {x(°\ x^- 1 ' , x^ 2 ' , ■ ■ ■ ,x^} of 
d + 1 d-cuts such that for each pair of j and k, the two d-cuts x^> and 
are adjacent. We call it an affine discrete cake-cut if we further require that 
Vj, k € {0, 1, • • ■ , d}: x^> and are affine adjacent. 

Definition 2. Discrete ENVY-FREE CAKE CUT: A discrete cake cut is 
an envy- free solution if there is a permutation ir of {0, 1, ■ ■ • , d} such that 
player i prefers the ir(i)-th segment for some d-cut x^' in the set. We denote 
it by Pi(xU)) = 7r(t). 

Note that, our definition is inspired and in line of the definition of discrete 
BROUWER fixed point [HI El [8]. 

Definition 3. 2D BROUWER: The input is a 2D grid of size G = N x N 
(N = 2 n ), together with a function f : G — > {0,1,2} such that a boundary 
condition is satisfied: \/y > : /(0, y) = 1, Vx > : f(x,0) = 2, and 
Vx,y > : f(x,N) = f(N,y) = 0. The required output is a unit square 
US = {(x, y), (x,y+l), (x + 1, y), (x + l,y+l)} such that f(US) = {0, 1, 2}. 

Using Kuhn's triangulation, labeling each node can be done in polyno- 
mial time, and so is coloring if the utility functions are given by a polynomial 
time algorithm. Therefore, each base cell can be constructed and their col- 
ors verified in polynomial time. In addition, vertices of each base cell in the 
Kuhn's triangulation are adjacent to each other. Therefore, the problem 
reduces to one of finding a fully colored Sperner cell, which can be done in 
PPAD. Therefore, we have the following: 

Corollary 1. Finding a discrete d-cut set for Envy-Free Cake Cutting 
problem for d + 1 people is in PPAD. 

On the other hand, we apply a reduction from the 2D BROUWER prob- 
lem to prove it PPAD-hard [7j; and hence: 

Theorem 3. Finding an approximate solution for Envy-Free Cake Cutting 
with d cuts for d + 1 people is PPAD -Complete. 

Proof. Given an input function of 2D BROUWER on grid / : N x N — > 
{0,1,2}, we embed it into a Kuhn's triangle defined by three vertices: 
< (0,0) J (2JV,0),(2JV,2JV) > by the mapping: M(x,y) = (2N - x,y). We 
define the preference functions (for all i = 0, 1, 2): Pi(2N — x, y) = f(x, y) for 
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< x,y < N; Pj(x,0) = 2, for < x < N; Pi(2N,y) = for N < y < 2N; 
P(x,y) = for all other cases. We name (0,0) the vertex X2, (2iV, 0) X±, 
and (2N,2N) Xq. Therefore, the Kuhn's triangle and the preference func- 
tions form a discrete ENVY-FREE CAKE CUT problem. The boundary 
condition for the SPERNER is now satisfied and there is a Sperner colored 
triangle, which is at the same time a ENVY-FREE CAKE CUT by our 
choices of the preference functions. Therefore ENVY-FREE CAKE CUT 
does have a solution. Once we find one, it must be in the region bounded 
by N < x < 2N, < y < N. The inverse mapping of M will give us the 
required BROUWER's solution. □ 

Therefore, the envy-free cake-cutting has the same time complexity as 
the Sperner Simplex computation if the utility functions are given by a 
polynomial time algorithm. 

Similar, under the oracle model for utility functions, we should show the 
same also holds. 

Theorem 4. Solving the ENVY-FREE CAKE CUT problem of d+1 people 
for the oracle functions requires time complexity 0{(^-) d ~ 1 )- 

Proof. By Lemma 2, the Kuhn's triangulation in the last section allows 
us to find a Sperner 's simplex in time 0{(^-) d ~ l ). Therefore, the solution 
corresponds to a discrete ENVY-FREE CAKE CUT solution. 

To prove the lower bound, we apply the same reduction for the Sperner's 
problem. Given an input instance of the zero point problem with direction 
preserving functions, we introduce the same structure as before. The only 
extra definition we need to introduce is for each vertex x of the grid, we 
define Vz : Pi{x) = j if f{x) = ej and Pi{x) = otherwise. Clearly, if 
a unit square contains all the preferences, there must be some i such that 
Pi(x) = 0. By direction preserving property, f(x) cannot be — ej for any j. 
Therefore f(x) = and a zero point is found. 

Our lower bound of follows. □ 

6 A Fully PTAS for 3 Players with Monotone Util- 
ity Functions 

For general oracle utility functions, the above results imply a 0(^-) matching 
bound when the number of players is three. Further improvement can only 
be possible when we have further restrictions on the utility functions. In 
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this section, we assume the utility functions in addition to satisfying the 
nonnegativity condition and the Lipschitz condition are also monotone. 

Even in this case, the celebrated Stromquist's moving knife [26] for envy- 
free cake cutting of three players can be shown to have an exponential lower 
bound (See Appendix A. 11). The main result here is to give an algorithm 
with a running time polynomial in the number of bits of K and -. 

Theorem 5. When the utility functions satisfy the above three conditions, 
an e envy- free solution can be found in time 0(log 2 when the number of 
players is three. 

The improvement has been made possible by an efficient way to compute 
the index of a triangulated polygon, i.e., counting the sum of the signs of 
base intervals on the boundary. The main idea is an observation that at 
some appropriate cuts, the colors along the cut are monotone. Therefore, 
we can find the boundary of the three different colors on this cut to calculate 
the sum quickly. 

At any time in the algorithm, we maintain a subset of V, with a non-zero 
index value, V(ii, 12, k%, fo) = {(hjj k) : i, j, k > 0,i + j + k = N,ii < i < 
i2,k± <k< k2}, delimited by i = ii,i = i2 and k = ki,k = k2- 

Algorithm 1. 1. If %i — i\ = 1 & k2 — k% = 1, find a fully colored base 
triangle ofV{i\,i\ -\-\,k\,k\ + 1), terminate. 

2. Choose max{i2 — i\, fc 2 — k\} (and assume it is 12 — i\ w.l.o.g.). 

3. Leti 3 = L(ii+* 2 )/2j. 

4- Calculate index (V(ii,«3, k\, A^)) and index(V(is, 12, k%, ^2))- 

5. Recurse on one of the two sub-polygons with a non-zero index. 

By definition, index(V(ii,t2, k\, ^2)) = index(V(ii,is, fei, k2))+index(V (i^, i2, k\ 
At least one of the index (V 13, k\, fe 2 )) and index{V{i^,i2, k±, fe 2 )) is non- 
zero as index(V(ii, 12, k%, ^2)) is non-zero. The algorithm always keep a 
polygon of non-zero index as the initial polygon has a non-zero index. At 
the base case ^2 = i\ + 1 and k2 = k\ + 1, V(i%, i\ + 1, ki, ki + 1) is either 
already a base triangle (with non-zero index) , or a diamond shape consisting 
of two base triangles, one of which must be of index non-zero. The correct- 
ness follows as the only possibility for its index being non-zero is when it is 
colored with all three colors. 

For complexity analysis, we first characterize the boundary conditions: 
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Property 2. There are up to three types of boundaries for V(ii,i2,k\,k2) 

1. Bi =c = {(i,j,k) € V : i constant, k\ <k< k2,i+j + k = N,i,j,k > 0} 

2. B k=c = k) € V : k constant, i\ < i < i 2 , i +j + k = N,i,j,k > 0} 

3. B j=0 = {(», 0, k) e V : i + k = N, i, k > 0} 

Second, we establish the monotonicity property. Note that for the third 
type of boundaries listed above, the monotonicity property does not hold 
for Bj =c in general but only for Bj = q. 

Property 3. The colors ofVoPiBi =c are monotone in k, and so are V\P\Bi =c 
and V2 H Bi =c . The same hold for R B k=c , as well as Vt H Bj = q, t = 0, 1, 2. 

Proof. Observe that, at Vo D -Bi= C! the color is determined by Player by 
finding the maximum of uq([0, ^]), uq([^, ^]), and Uo([^, 1]), where < 
j < N — c. The first item is fixed. The second item is increasing in j 
since [fj, C [-^, ^r~] for j < / and m is assumed to be a probability 
distribution. For the same reason, the last item is decreasing in j. The color 
will be if «o([0, jj]) is the maximum, 1 if Uo([jt, ^ff-]) is the maximum, or 
2 otherwise. In general, as j increases, the color in Vo H Bi =c will start with 
2, then 0, and finally 1, assuming non-degeneracy. However, any of those 
colors may be missing. 

The same analysis holds for other players and for the case when k is 
fixed. However, it does not hold when j is fixed, except when j = 0. □ 

Property[3]allows us to find the colors of all the vertices on the boundaries 
V t riB i=il ,V t nB i=i2 , V t nB j=0 , V t nB k=kl and V t nB k = k2 , t = 0, 1,2, in time 
proportional to logarithm of the length of the sides, by finding all vertices 
at which the color changes along it. Once the changing points of the colors 
are found, using monotonicity, the total number of positively and negatively 
signed base intervals along the boundary can be calculated in constant time. 

The recursive algorithm reduces the size of V(ii, i%, ki, A^) geometrically. 
Let L = max{|ii — za | , \ki — foQ- L halves in two rounds of the algorithm. 
It takes 21og 2 -/V steps to reduce L to 1. 

At each round of the algorithm, we need to find out, for each of three 
players, for each side of V(ii, 12, k\, A^) (up to five in all), the boundary of the 
player's color changes (two boundaries). That takes time log 2 L, bounded 
by log 2 N, to derive a total of 3 x 5 x 2 log 2 N. 

Therefore, the time complexity is 0((log 2 N) 2 ). 
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7 Discussion and Conclusion 



It remains open whether the approximate envy-free cake cutting problem 
for four or more players would allow for a polynomial time approximation 
scheme as in the three player case, when we are dealing with measurable or 
monotone functions. 
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Appendix 

A.l Figure [D Base triangle with edge orientation 



(3.0.0) 




(0,3 : 0) (0.2,1) (0,1.2) (0,0,3) 



Figure 1: Base triangle with edge orientation 
A. 2 Proof of Lemma Q] 

Proof. By Proposition [1] and the definition of index(A, 4>), the statement 
that there are at least index(A, 0) fully colored base triangles is obviously 
true. 

Note that every internal base edge is in two base triangles and has dif- 
ferent signs with respect to the two base triangles. Therefore, they cancel 
each other and to derive the following: 

index(A, (ft) = ^^{sign(e, (ft) : e a boundary base edge} 
Therefore, we complete the proof. 

Note that in the above sum, the orientation of the boundary base edges 
are in the clockwise order around A. □ 

A. 3 Proof of Proposition [2] 

Proof. For each cf-dimensional simplex index<i(5, (ft) is one if and only i f the 
set of colors of its d+1 vertices are all distinct and is the same as {0, 1, ■ • • , d}. 
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On the other hand, it has d+ 1 faces of d — 1-dimension. Each face of (d— 1)- 
dimension is a simplex of (d— l)-dimension. By the definition of the induced 
index indexd-i only considers colors {0, 1, • • • , d — 1}, a [d — l)-dimension 
simplex has an index 1 if and only if the set of colors of all its vertices is 
the same as {0, 1, • • • , d— 1}. Now there is only one vertex left for which we 
don't know its color. If it is anything in {0, 1, • • • ,d— 1}, we have exactly 
one more face having index equal to 1. In that case, summing up the indices 
of the faces, we obtain a sum of two which is zero mod 2. If the last color is 
d, then the sum will be one. Therefore, the claim holds when P is a simplex 
with no further triangulation. 

In general, index(P,(p) = YlseP indexes, <f>). We can replace index(8, (ft) 
by the sum of indices of the boundaries of 5. However, each (d — 1)- 
dimensional simplex in the triangulation appears in exact two d-dimensional 
simplices in the triangulations, unless it is at the boundary of P where it 
appear only once. Since we consider the sum mod 2, all the terms cancel 
out except those on the boundaries of P. The claim follows. □ 

A. 4 An example of Kuhn's triangulation for a 3-cube 

Let (0, 0, 0) be the base point, according to the different permutations of 
0,1,2, we obtain six tetrahedrons which is a simplicial partition of the unit 
cube. 



7T = 


(0,1,2) 


A = 


{(0,0,0), (1,0,0), 


(1,1,0), 


(1.1,1)} 


7T = 


(0,2,1) 


A = 


{(0,0,0), (1,0,0), 


(1,0,1), 


(1,1,1)} 


7T = 


(1,0,2) 


A = 


{(0,0,0), (0,1,0), 


(1,1,0), 


(1,1,1)} 


7T = 


(1,2,0) 


A = 


{(0,0,0), (0,1,0), 


(0,1,1), 


(1,1,1)} 


7T = 


(2,0,1) 


A = 


{(0,0,0), (0,0,1), 


(1,0,1), 


(1,1,1)} 


7T = 


(2,1,0) 


A = 


{(0,0,0), (0,0,1), 


(0,1,1), 


(1,1,1)} 



The partition is illustrated in Figure 
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Figure 2: An illustration of Kuhn's partition on a unit cube in 3 dimension 



A. 5 Proof of the consistency of the two processes 

we scale the unit cube into one of side length N. We call it the big cube, 
and each of the N d sub-cubes as the small cubes. 

It is enough to show that, for a simplex in the refined grid starting from 
any grid point x* = (a^a^, ■ ■ ■ , x£) to x* = (x\ + 1, x\ + 1, ■ • • , x* d + 1), 
defined by a permutation 7, all its d+1 vertices are in the same big simplex 
derived by some permutation p with base point (0, 0, • • • , 0) to the point 
(N,N,--- ,7V). 

Let x° = x*, and x l is the same as x l ~ l except in its coordinate x^^ 
which is x l ~,\ + l. Therefore, if x* > x*, then x* > x* + l and then x k > x k for 
all k = 1, 2, ■ • ■ ,d. Similarly, if x* < x*, then x\ < x 1 * for all k = 1, 2, • ■ ■ , d. 
If x* = Xj, then we have x\ > x k - for all k = 1, 2, • ■ ■ , n, or x\ < Xj for all 
k = 1,2, • • • , n, dependent on 7(1) < 7(j) or 7(1) > 

Therefore, there is a permutation p such that x k p ^ > x k p ^ >••• > x pid) 
for all k = 0, 1, 2, • • • , d, which guarantees the base simplex inside one of the 
d\ large simplices. 



A. 6 Transformation process 

This can be done by a transformation as follows. Let X be a base point 
in barycentric coordinates. Set = (eio, en, ■ ■ ■ , e^) where = for all 
i and j except that en = —1, e^j+i = 1 for z = 0, ■ • • , d — 1. Then the 
vertices of the base simplex according to permutation tt based on X are 
given by v\ = v % ~ 1 + e^ (j) for Vz = 1, • ■ ■ , d and = X. For example, for 
a simplex based on vertex v = (0, 0, • ■ ■ , 0)^ corresponding to permutation 
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(0, 1,2, ••• ,d— 1). We first set v° = (1,0, ••■ ,0)d+i, then by the above 
transformation v 1 = (0, 1, • • • , 0)d+i, v 2 = (0, 0, 1, • • • , 0)^+1- Please refer to 
page 42 of [23j for more details. 

A. 7 Proof of Lemma [5] 

Proof. WLOG, assume the base cell corresponds to a permutation tt and 
Y = X + J2i=i e 7Tj f° r some < I < k < d. Fix i, since en = 1, e^j+i = — 1 
are the only non-zero coordinates in ej, let a = (oi, a-2, • • • , = X^=z e 7r;> 

then Vz, aj £ {—1,0, 1}. Therefore, 5x-y equals to either or 1 and since 
X, Y are two different vertices, we must have 5x-y = 1- D 

A. 8 Proof of Lemma [4] 

Proof. We establish a reduction from the zero point problem considered by 
Chen and Deng [6] for the direction preserving functions, to the Sperner 
problem of dimension d to derive a similar matching bound within a constant 
factor for any constant dimension d. 

We achieve the goal in two steps. Let N = |. First, we reduce the 
problem in [6] to one of finding a discrete fixed point over a base hypercube 
on hypergrids for functions from the grid points N d to {0, 1, 2, • • • , d} sat- 
isfying the usual boundary conditions, where a discrete fixed point over a 
base hypercube (2 d points within distance one m J • |qq metric) is one with 
function values on the nodes of the base hypercube cover all the values from 
to d. 

To do that, we quickly review the fixed point problem for direction 
preserving functions discussed in [6]. We consider functions which are de- 
fined on the grids point N d , which has values Ud := {0, ±ei, ±e2, ■ • • , ie^}, 
where is the unit vector with the i-th. coordinate being 1. A function 
/ : N d — > Ud is direction preserving if and only if for each x, y S N d with 
1^ — y\oo < 1> f(x) T ■ f(y) > 0. That is, the function values of nodes within 
distance 1 of each other cannot be of different signs. We define a function 
h : N d -> {0, l,--- ,d} such that h(x) = if f{x) < and h(x) = i if 
f(x) = ej. If we find a discrete fixed point set in h on a base hypercube, 
then the base hypercube must contain a node x such that f(x) = as / 
is direction preserving. Therefore, finding a zero point in / can be reduced 
to the problem of finding a discrete fixed point of the hypercube form in 
function values of h. There is still an issue whether such a solution exists. 
For the zero point problem, in addition to the requirement of direction pre- 
serving, it is assumed that the function / is bounded, i.e., f{x) + x G N d . 
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In fact, we can further assume a specific boundary condition that either 
f(x) = if Xi = and Vj < i : Xj > 0, f(x) = —ei if x » (all the 
coordinates of x are positive), Xi = N and Vj < i : Xj < N. In fact, for 
any bounded function /, we can always cover it with one more layer on 
each face of its boundary (and function values) to achieve that. It is easy 
to verify the direction preserving condition still holds for bounded functions 
satisfying the property in the original cube. Now it is easy to see there is 
only one cube of dimension {d — 1) that has an index one on the boundary. 
By Proposition [2j there must be a base simplex inside the hypercube that 
has index one, i.e., has its vertices colored differently from {0, 1, • • • ,d}. 



(0.2) (1.2) [2.2] 




Figure 3: Reduction process in Lemma 3 

In the second (and last) step, we note that there is a size (^r) d hypergrid 
in a d dimensional simplex of length N. The lower bound ^((-j)^ 1 ) follows 
from a lower bound of Q(N d ~ l ) for the fixed point problem. Since d is a 
constant here, we obtain the lower bound f2(iV d_1 ). We need to add that, 
the hypergrid can be embedded into the simplex such that its coloring can 
be extended to a valid coloring of vertices in the triangulated simplex. We 
align the colored hypergrid in a way that the origin is placed at the origin 
of the simplex, and align the d rays out of the hypergrid with the rays of 
the simplex out of its origin. For the overlapping vertices of the simplex, 
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they are colored with the same colors as in the hypergrid. The faces of the 
simplex passing through the origin will be colored by the same rules as the 
corresponding faces of the hypergrid, except the vertices Ne\, Ne2, • ■ ■ , Ned- 
which are colored in the following rules: Nei is colored with i + 1 and Ned 
is colored with 0. Note that the origin is colored with 1. All the extra 
vertices of the simplex will be colored 0. This way, the boundary conditions 
of the simplex are satisfied. In addition, no new fully colored base simplex 
is introduced into the triangulated simplex. See Figure [3l 

□ 

A. 9 Figure [4} Barycentric subdivision on a two dimensional 
triangle 




(a) 1st iteration (b) 2nd iteration (c) 3rd iteration 



Figure 4: Barycentric subdivision on a two dimensional triangle 
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A. 10 Figure [5} An illustration of Sperner simplex approach 
for 3 players envy-free cake-cutting 




Figure 5: Sperner simplex approach for 3 players envy-free cake-cutting 

A. 11 Proof of an exponential lower bound for Stromquist's 
solution 

The celebrated Stromquist's solution [26J involves a referee who moves her 
sword from left to right. The three players each has a knife at the point 
that would cut the right piece to the sword in half, according to their own 
valuation. While the referee's sword moves right, the three knives all move 
right in parallel but possibly at different speeds. At all times, each player 
evaluates the piece to the left of the sword, and the two pieces that would 
result if the middle knife cuts. If any of them sees the left piece of the sword 
is the largest, he would shout "cut". Then the sword cuts. The leftmost 
piece is assigned to the player who shouted. For the two players who didn't 
shout, one whose knife is to the left of the middle knife receives the middle 
piece, and one whose knife is on the right of the middle knife receives the 
rightmost piece. 

we will show that Stromquist's moving knife procedure can not be turned 
into a polynomial time algorithm for finding an e-envy free solution. We will 
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do this by showing that the particular fixed point found by the Stromquist's 
procedure can not be found with polynomial number of queries. We assume 
a query can ask about the valuation of a player for an interval (x,y). 

Suppose we have three players A, B, and C. The cake is represented as 
an interval [0,1]. Players A and B have the same utility function. Their 
utility functions can be described as follows q 



ua{x,v) = u B (x,y) 



for x = 0,y = 1/10 

2 for x = 1/10, y = 3/10 

100 for x = 3/10, y = 4/10 

2 for x = 4/10, y = 8/10 1 

100 for x = 8/10, y = 9/10 

for x = 9/10, y = 1 

The value of both players for any other interval can be computed as- 
suming that their valuation is uniform across all the above intervals. For 
example (1/20, 3/20) = 0.5. uc can be described in a similar way as 
follows: 

100 - 5 for x = 0,y = 1/10 
2 for x = 1/10, y = 3/10 

uc(x,y) = <j 98 
2 




_ ^ ? — / 

for x = 3/10, y = 4/10 
for x = 4/10, y = 5/10 
for x = 5/10, y = 1 



Assume that the value of 5 is very small. Now, consider the Stromquist's 
moving knives. Suppose the referee starts by moving her knife from to- 
wards 1. A and B have the same utility function so their knives are going 
to be at the same place. Moreover, the "middle knife" will be always A's 
or B's. Observe that when the referee's knife reaches point 1/10 the middle 
knife is at point 4/10. No one will shout cut yet. Now, observe that since 
according to A and B the density of the interval (1/10, 3/10) is twice the 
interval (4/10, 5/10), the speed of the knives of A and B will be exactly the 
same as the referee's knife. For a similar reason, according to C the value of 
the middle piece (from referee's knife to the knife of A or B) will remain just 
5 above the value of the leftmost piece until the referee's knife goes slightly 
beyond 3/10. At that time C will shout cut and receive the leftmost piece. 
Players A and B happily split the rest of the cake equally. 



1 We should normalize them for consistency but did not for simplicity of presentation. 
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Now, we will slightly perturb the utility function of C. At a point 1/10 < 
x < 3/10 perturb the utility function of C in the following way: increase the 
utility of C for the interval (x — 5, x) by 5/2 and decrease his utility for the 
interval (x, x + <5) by the same amount. Call the player with this perturbed 
utility function C x . 

Now, it is easy to see that if we are running the Stromquist's method, 
player C x will shout cut at the time referee's knife reaches x. It is also 
not hard to see that any value query asking for the value of an interval can 
not distinguish C from C x unless one end points of the interval for which 
it is querying the value is in {x — 8, x + <5). Therefore, it is not possible to 
distinguish C x from C with queries of order polylogarithmic in 1/8. 
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